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ABSTRACT 

A perturbation  formulation  of  the  equations  of  linear  piezoelectricity 
for  small  fields  superposed  on  a bias  is  obtained  from  a Green*  s function 
representation.  It  is  shown  that  the  resulting  equation  for  the  first  per- 
turbation of  the  eigenvalue  may  be  obtained  without  the  use  of  a Green*  s 
tensor  or  a complete  set  of  orthogonal  eigensolutions . Since  the  bias  enters 
the  constitutive  equations,  the  boundary  conditions  contain  perturbation  terms 
as  well  as  the  differential  equations.  The  linear  electroelastic  equations 
for  small  fields  superposed  on  a bias  differ  from  the  equations  of  linear 
piezoelectricity  because  the  effective  material  constants  of  linear  electro- 
elasticity have  less  symmetry  than  the  constants  of  linear  piezoelectricity. 
Consequently,  a perturbation  formulation  of  the  linear  electroelastic  equa- 
tions for  small  fields  superposed  on  a bias  is  presented,  it  is  shown  that 
the  effective  constants  of  linear  electroelasticity  have  just  the  symmetry 
required  for  the  condition  of  the  orthogonality  of  linear  electroelastic 
vibrations  to  hold. 
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1.  introduction 

In  vibrating  piezoelectric  solids  small  chemges  in  such  things  as  natural 
frequency  emd  small  values  of  attenuation  aurising  from  a variety  of  causes  may 
readily  and  conveniently  be  determined  from  a perturbation  formulation^  ^ of  the 
equations  of  linear  piezoelectricity  by  employing  the  equation  for  the  first 
perturbation  of  the  eigenvalue,  liiis  fact  is  well  known  and  widely  appreciated. 
The  small  effects  may  be  due  to  such  things  as  material  viscosity,  air  loading 
of  the  surface,  mass  loading  and  stiffening  of  the  surface  due  to  the  presence 
of  a thin  surface  film,  the  electrical  conductivity  of  the  film  and  biasing 
stresses,  strains  amd  electric  fields,  which  either  may  be  applied  or  a conse- 
quence of  such  things  as  changes  in  temperature.  Although  the  linear  piezoelec- 
tric equations  may  readily  be  employed  in  the  treatment  of  most  of  the  above- 
. mentioned  phenomena,  they  are  invalid  when  the  above-mentioned  biasing  states 

are  present  because  the  linear  piezoelectric  equations  are  not  equivalent  to 

* 4 

the  linear  electroelastic  equations  for  small  fields  superposed  on  a bias  . The 

basic  reason  for  this  is  that  although  the  linear  piezoelectric  equations  are 

the  appropriate  linear  limit  of  the  properly  invariant  nonlinear  electroelastic 

equations^,  they  are  not  the  appropriate  linear  limit  of  the  properly  invariant 

nonlinear  electroelastic  equations  for  small  fields  superposed  on  a bias. 

In  this  paper  perturbations  from  solutions  of  the  linear  piezoelectric 

equations  due  to  the  above-mentioned  biasing  states  are  treated  by  means  of  the 

aforementioned  linear  electroelastic  equations  for  small  fields  superposed  on  a 

bias.  In  the  treatment  the  equations  are  written  as  the  linear  piezoelectric 

equations  plus  the  additional  terms  arising  from  the  bias,  and  a piezoelectric 

perturbation  theory  is  obtained  from  a Green' s function  formulation  of  the 

equations  of  linear  piezoelectricity,  ihe  derivation  is  identical  with  that 

2 

given  in  a recent  work  , but  the  perturbation  terms  are  different  and  only 


y 
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perturbations  due  to  biasing  states  are  considered,  ■flie  small  effects  due 
to  the  aforementioned  other  phenomena  are  expressly  not  considered  here,  but 
the  inclusion  of  them  would  be  straightforward  and  may  readily  be  accon^lished 
after  the  fact.  The  resulting  equation  for  the  first  perturbation  of  the 
eigenvalue,  which  is  also  obtained  in  a completely  independent  manner  withouc 
assuming  the  existence  of  a complete  set  of  orthogonal  eigensolutions,  may 
readily  be  applied  in  the  determination  of  small  changes  in  natural  frequen- 
cies of  vibrating  piezoelectric  solids  subject  to  various  biasing  states. 
Indeed,  the  perturbation  equation  has  already  been  applied^  in  the  calculation 
of  the  change  in  velocity  of  piezoelectric  surface  waves  on  Y-cut,  Z-propa- 
gating  lithium  niobate  subject  to  uniaxial  biasing  stresses.  The  results  of 
the  perturbation  calculation  agree  so  well  with  complete  calculations  performed 
for  the  same  case  that  it  may  be  concluded  that  the  perturbation  procedure  is 
every  bit  as  accurate  as  the  complete  calculation  for  the  determination  of 
such  quantities.  In  fact,  since  the  perturbation  procedure  can  readily  treat 
spatially  varying  biasing  states,  for  which  a complete  calculation  cannot  be 
performed,  it  has  a significant  advantage  over  a complete  calculation  in  the 
determination  of  this  type  of  quantity. 

Since  it  is  possible  for  perturbations  to  arise  from  the  above-mentioned 
other  phenomena  even  when  a biasing  state  is  present  and  the  piezoelectric 
perturbation  theory  is  invalid  under  such  circumstances,  an  electroelastic 
perturbation  theory  is  obtained  from  a Green' s function  formulation  of  the 
linear  electroelastic  equations  for  small  fields  superposed  on  a bias.  This 
treatment  differs  from  the  previous  one  in  that  the  terms  arising  from  the  bias 
are  included  in  the  basic  equations  that  must  be  satisfied,  rather  than  as  a 
perturbation  from  the  linear  piezoelectric  equations.  In  this  linear  electro- 
elastic case  the  effective  material  constants  have  less  symmetry  than  in  the 


3. 


linear  piezoelectric  case,  but  they  have  just  the  synrietry  required  for  the 
orthogonality  theorem  for  linear  electroelastic  vibrations  superposed  on  a bieis, 
which  we  need  for  our  perturbation  theory,  and,  of  course,  the  related  theorems 
of  reciprocity  and  vmiqueness.  The  resulting  equation  for  the  first  perturba- 
tion of  the  eigenvalue  may  readily  be  applied  in  the  determination  of  small 
changes  in  natural  frequency  and  phase  velocity  and  small  values  of  attenuation 
of  vibrating  biased  electroelastic  solutions,  which  arise  from  what  we  have 
termed  other  phenomena. 

2 . Linear  Electroelastic  Equations  for  Small  Fields  Superposed  on  a Bias 

Before  presenting  the  equations  we  briefly  introduce  some  preliminary 

terminology  and  notation.  We  first  note  that  under  the  static  bias  the  material 

points  move  from  the  reference  coordinates  X,  to  the  intermediate  coordinates  ? 

L 

and  we  have 

= (2.1) 

Then  in  the  superposed  small  dynamic  motion  the  material  points  move  from  the 
intermediate  coordinates  to  the  present  coordinates  y^,  and  we  have 

^i  " • (2.2) 

We  consistently  use  the  convention  that  capital  Latin  indices,  lower  case 
Greek  indices  and  lower  case  Latin  indices,  respectively,  refer  to  the 
Cartesian  components  of  the  reference  coordinates,  intermediate  coordinates 
and  present  coordinates  of  material  points.  In  this  paper  cartesian  tensor 
notation  is  used  exclusively.  A coitma  followed  by  an  index  denotes  partial 
differentiation  with  respect  to  a geometric  coordinate,  i.e.,’ 
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'c.L  ■ i'i,.  - - 

%i  - h,t  ■ \,i  ' . 


(2.3) 


and  we  employ  the  summation  convention  for  repeated  tensor  indices  and  the  dot 
notation  for  differentiation  with  respect  to  time,  since  the  dynamic  motion 
is  small,  we  may  write 

yi  = ^ip<%  + , (2.4) 


where  u-  is  the  small  mechanical  displacement  from  the  intermediate  coordinate 
P 

to  the  present  coordinate  y . and  6 . „ is  a Kronecker  delta,  which  serves  to 

X xp  ^ 

translate  a vector  from  the  present  to  the  intermediate  coordinates  and  vice 
versa. 

By  virtue  of  (2.1)  the  small  field  dynamic  equations  may  use  either  the 
intermediate  coordinates  5^  or  the  reference  coordinates  x as  independent 
variables.  Since  in  the  typical  situation  it  is  undesirable  to  measure  the 
geometry  each  time  the  bias  is  varied,  it  is  advantageous  to  use  the  X as 
independent  variables.  Moreover,  in  the  applications  envisaged  although  the 
static  biasing  mechanical  displacement  w of  a material  point  is  large  compared 
with  Up,  it  is  nevertheless  small,  and  we  may  write 


•'a 


(2.5) 


where  6 is  a Kronecker  delta,  which  serves  to  translate  a vector  from  the 
intermediate  to  the  refezence  coordinates  and  vice  versa  and  is  required  for 
notational  consistency  amd  clarity  because  of  the  use  of  different  indices 
to  refer  to  the  different  coordinates  of  a material  point.  In  accordance  with 
the  foregoing  the  present  electric  potential  at  a material  point  may  be 
written 


^ 1 

n 


cp(Yj,t)  = 9(Xj^,t)  = 9 (Xj^)  + 9(Xj^,t)  , 


(2.6) 


where  cp^  is  the  biasing  electric  potential  at  a material  point  and  9 is  the 
small  field  dynamic  electric  potential  at  the  same  point. 

since  we  have  used  the  reference  coordinates  X as  the  independent  variables 

Xj 

for  material  points,  we  take  the  X^^  rather  than  the  as  the  independent 
variables  for  points  of  free  space  not  abutting  a material  boundary.  Moreover, 
since  only  the  reference  position  of  the  material  boundary  is  known,  the  value 
of  the  electric  potential  t on  the  free  space  side  of  the  present  position  of 
the  material  boundary  must  be  obtained  by  means  of  a Taylor  expansion  about 

its  value  at  the  known  reference  position  of  the  material  boundary  in  the 
unknown  mechamical  displacement  (w^  + of  material  boundary.  Thus 

nyyt)  = l^x^)  t ^^^W^  + ^ 

''1  . . . 

to  first  order  in  the  small  field  variables,  cuid  where  lit  is  the  biasing 
electric  potential  immediately  on  the  free  space  side  of  the  reference  position 
of  the  material  boundary  and  \|t  is  the  small  field  dynamic  electric  potential 
at  the  same  point  and  for  points  of  free  space  we  have 


t (X  t)  = t (X  ) + t(X  ,t)  . 


(2.8) 


Now  that  the  meaning  of  the  basic  dependent  and  independent  variables 
has  been  explained  we  are  in  a position  to  write  the  linear  electroelastic 

4 

equations  for  small  fields  superposed  on  a bias,  which  take  the  form 


where 


iL'NMo'^ Qf,M 
\ Of,M 


(2.10) 


ILVMD'  ^L’VMa  2L’VMQ''  ful  ” ®LM  ” ®LM' 


G " R ~ G ■f'  © 

2ML-Y  iML-y  MLY  MLY  ’ 

- 1 . 1 
c =t6+c  E+c  w 

2LVMa  LM  'la  3LYMa^AB  AB  2li'YKM  «,  K 

, -1  1 

+ c w -k  cp+g 

2LKMa  Y,  K lALYMa^^A  ^LYMcy, 

- 1 a1  2 

e =-k  E-ew  -b  cp+g  , 

MLY  IMLYBC  BC  MLK  Y,  K AMLY  ,A  MLY  ’ 

1 - 1 ^ 1 J. 

e =b  E -X_cp  -2ejE. 

LM  LMCD  CD  ^MC  % C O ML’ 


(2.11) 


(2.12) 


and 


A,  Ac 

" 4®"a,  A,v  - ' 

^MLY  " ^ ^o^^0\p\,Y  " VL,a\Y  " ^Y^M,c^\,cy^  ’ 


(2.13) 


and  we  have  taken  the  liberty  of  changing  the  order  of  the  last  two  indices  on 

the  second  and  third  rank  tensors  having  both  Latin  and  Greek  indices  from  that 

employed  in  Ref. 4.  Ecjuations  (2.9)  are  the  small  field  stress  equations  of 

motion  and  charge  ecjuation  of  electrostatics  referred  to  the  reference  ccjordi- 

nates  X . Equations  (2.10)  are  the  linear  electroelastic  constitutive  equations 
L 

and  Eqs. (2.11)  and  (2.12)  contain  the  definitions  of  the  effective  coefficients 

'ft 

defined  therein.  in  (2.9)  - (2.13)  and  denote  the  components  of  the 

small  field  Piola-Kirchhoff  stress  tensor,  which  is  asymmetric,  and  reference 

4 

electric  displacement  vector  , respectively;  denotes  the  reference  mass 

density,  c ..  . e and  e denote  the  second  order  elastic,  piezoelectric 
2LYMa’  MLY  LM  ’ ^ 

and  dielectric  constants,  respectively,  which  are  the  constcintp  that  occur  in 

the  ordinary  linear  theory  of  piezoelectricity.  The  symbols  T^,.  and  E^_  denote 

LM  AB 

the  conponents  of  the  static  biasing  stress  and  strain,  respectively  and 
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is  the  Jacobian  of  the  static  deformation.  Hie  biasing  variables  satisfy  the 
appropriate  static  equations  given  in  (66)  - (72)  of  Ref. 4,  or  the  equivalent 
equations  using  reference  coordinates  as  independent  variables.  Since  we  are 
interested  in  small  biasing  strains  only,  we  have 


4 (w  + w ) . J^  ^ 1 
AB  2 ^ A. B B.A^  ’ ’ 


1 

T = c w + e CP 
LM  2LMRS  R,  S RLlT,  M ’ 


(2.14) 

(2.15) 


and  each  X in  (2.13)  may  be  replaced  by  6 . In  (2.12)  c . b 

1^,01  y i LQ;  ' ' 3LYMQAB’  AMLV^ 

and  denote  the  third-order  elastic,  electros  <.rictive,  third-order 

electric  permeability  and  first-order  electroelastic  co’-  stants,  respectively, 
and  denotes  the  electric  permittivity  of  free  space.  For  obvious  reasons 
the  notation  employed  here  is  designed  to  be  consistent  with  the  notation  of 
Ref. 4.  The  karets  over  many  variables  have  been  employed  here  because  we  con- 
sistently use  the  reference  coordinates  as  independent  variables. 

To  the  foregoing  ecjuations  we  must  adjoin  the  dynamic  boundary  conditions, 
which  when  the  boundary  of  the  body  abuts  free  space  take  the  form 

? ‘ * *,1,^  ^ ’ 


(2.16) 


where  w^  and  u^  in  (2.16)2,  denote  the  static  biasing  and  small 

field  dynamic  components  of  the  mechanical  displacement  at  the  surface  of  the 

solid,  N denotes  the  components  of  the  unit  normal  to  the  reference  position 

of  the  surface  and  T^  is  an  applied  traction  per  unit  reference  area.  Ihe 

free-space  dynamic  variables  K cind  3 are  given  by 

liV  L 
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'''  2LVMQ;'^a,M  fLyM*,M  ’ 


2^  = 


L = jLg'^3  ^ ?LmAm  ^ ? ^M  ' 


(2.17) 


where 


fHS  ■ * '^vya  - ’ 


fLv«c“  A,  b'o‘-=^Xc  * I 


\»%>  * =J?Vcb'  ’ 


Imm*  A,b'o<=^«V  - ''v‘«B  - =i'vB>  > 


f 1 

= J X. 


5mv  ■ ■’A,  A,  B'o'VvB  - '^b'V  ’ 


f 1 1 „f  1 

Rq  = JX  se  q,  R=-Je  , 
£lP  L, Of  o O'  P’  3 o 


(2.18) 


In  free-space  the  small  field  dynamic  electric  potential  clearly  satisfies 


Laplace's  equation  ^ i.e., 


LL 


= 0 . 


(2.19) 


When  the  body  abuts  another  solid  insulator  rather  than  free-space,  instead  of 
(2,16)  we  must  have 

(2.20) 


If,  on  the  other  hand,  the  body  abuts  a perfect  conductor,  we  have 

9=^,  (2.21) 

~ ~ 

where  jS'^  in  the  conductor  vanishes,  in  the  conductor  is  not  given  by  (2.10), 


is  purely  mechanical  and  depends  on  the  particular  case  treated  and  Si  is  the 
dynamic  electric  surface  charge  density  per  unit  reference  area.  In  the  usual 
case  of  interest  cp  will  be  prescribed  and  (2.21)^  will  determine  the  surface 
charge  density  a posteriori. 


ij 
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In  this  section  we  obtain  the  equations  for  perturbations  from  solutions  of 

the  linear  piezoelectric  equations  due  to  a bias  from  a Green' s function  formu- 
2 

lation  of  the  equations  of  linear  piezoelectricity.  To  this  end  we  first 
write  Eqs.  (2.9)  - (2.12)  in  the  form 

= P°u  , (3.1) 

Ly,L  L\,L  V’ 

" 2L\MQ''^ OjH 

= e u - e CD  , (3.4) 

L LMOf  a,M  LM%M’ 


(3.3) 


(3.4) 


where  the  nonlinear  terms  ^ and  ^ are  the  perturbation  terms  which  depend 

L 1 Xj 

-'1 

on  the  biasing  state  w and  cp  in  accordance  with  (2.12)  and  (2.14)  . 

6 ^ M X 

since  K?  and  are  the  perturbation  terms  we  are  perturbing  frcjm  the  linear 
liT  Ij 

2 

piezoelectric  equations  and,  hence,  the  equations  for  the  Green' s tensor  may 


be  written  in  the  form 


K^^,L  tpcD  g“=-  6(P-Q)6“, 
^L,L=- 

~ia  a .Of 

^•Y  2LVMC®C,M 

= e ,g“  - e f“  , 

L LMC  C,M  LM  ,M’ 


(3.5) 


(3.7) 


where  the  superscripts  Qf,  P = 1 - 4,  P and  Q denote  the  fixed  field  point  and 

• Q 

varieible  source  point,  respectively,  6 is  the  Dirac  delta  function,  6^  is  the 
a o 

Kronecker  delta,  G^  and  f are  the  mechanical  displacement  Green' s tensor  and 
electric  potential  Green's  function  plus  cross  terms,  respectively.  In 


10, 


t I 

I ; 


i 

1 

i 

I 

! 

I 

I 


iu)t 

(3.5)  - (3.7)  we  have  assumed  that  all  variables  have  a time  dependence  of  e , 
which  has  been  factored  out,  emd  we  incOce  the  same  assumptions  in  (3.1)-  (3.4). 
Equations  (3.1)-  (3.3)  and  (3.5)-  (3.7)  are  identical  with  Eqs.  (2.1)-  (2.3)  amd 
(2.4)-  (2.6)  of  Ref. 2,  except  for  notation.  Consequently,  except  for  notation. 
Sec. 2 of  Ref. 2 applies  here  without  change.  Thus,  by  following  a procedure 
identical  with  the  one  employed  in  Sec. 2 of  Ref. 2,  from  Eq.  (2.16)  of  Ref . 2 with 
the  change  in  notation  we  may  write 


u^  = E -V ^ 

0^  II  ...2  _.2 


JJL  O 


[J  \ £ (2)  lel  - K (5)1  (8) 


(3.8) 


where  U)  is  the  eigenfrequency  of  the  p-th  eigensolution  and 


Y (M-)  ’ ^ ^“(p.)  ' 


(3.9) 


in  which  u[^  and  cp  denote  the  pth  eigensolution  functions  satisfying  (3.1)  and 

(3.2)  in  the  absence  of  K?  and  ^ and  siibject  to  the  appropriate  homogeneous 

Xj  V L 

7 

boundary  conditions  and  N , may  be  found  from 

(P ) 


i p u^  u''^  dv  = 6 

J ‘'o  Y Y o “(U)  Pv^ 


(3.10) 


; 

1 


■i 


which  is  the  orthogonality  condition  for  piezoelectric  vibrations.  iSius,  g*|^ 

and  f denote  orthonormal  eigensolutions  of  the  appropriate  homogeneous  problem. 

in  (3.8)  S denotes  the  reference  surface  area  enclosing  the  reference  volume  V 
o o 


and 


11. 


we  may  now  obtain  the  perturbation  procedure  from  (3.8)  in  the  usual  way  , 
i.e.,  by  letting  u^(P)  be  very  near  one  of  the  say  g^ . Ihen  we  may  write 


where 


" I - "A  * 

*1  ^'’o- 


(3.12) 


(3.13) 


and  from  (3.8),  (3.12)  and  (3.13)  we  have 

» 1,  (3.14) 

which  is  the  equation  for  the  first  perturbation  in  eigenfrequency.  if 

A = 1^1  (3.15) 

from  (3.14)  we  have 

A * Hj^2u)^ , (3. 16  ) 

for  the  first  perturbation  in  eigenfrequency.  We  have  discussed  first-order 
perturbation  theory  only  because  that  is  all  we  are  interested  in  here.  For  a 
discussion  of  second-auid  higher-order  perturbation  theory  see  Ref. 8. 

The  equation  for  the  first  perturbation  of  the  eigenvalue,  i.e.,  Eq.  (3.16), 
is  a very  importemt  relation  that  has  numerous  applications.  Although  it  was 
obtained  from  a Green's  function  formulation  by  assuming  the  existence  of  a 
con^lete  set  of  orthogonal  eigensolutions,  this  particular  relation  cem  readily 
be  obtained  without  the  use  of  a Green' s tensor  or,  more  importamtly,  a complete 
set  of  orthogonal  eigensolutions.  in  fact,  the  existence  of  only  the  particular 
unperturbed  eigensolution  under  consideration  and  a nearby  perturbed  state  is 
required.  To  see  this  consider  the  unperturbed  Mth  eigensolution,  at  eigen- 
frequency which,  from  (3.1)  and  (3.2),  satisfies 


I 


i 

! 


! 

i 

4 

I 

■ 

I 

! 


i(f“  *pV»"  = 0,  5*"  .0, 

LY,L  ^ M Y L,L 


along  with  the  nearby  perturbed  solution  at  frequency  ou,  which,  from  (3.1) 


(3.17) 


and  (3.2),  satisfies 


where  it  is  understood  that  (3.17)  and  (3.18)  eure  independent  of  time.  From 


(3.17)j^  and  (3.18)^^  form  the  equation 


Performing  the  usual  operations  , employing  (3.3),  (3.17)2  and  (3.18)2  and  the 
divergence  theorem  in  the  usual  manner,  we  obtain 

V s 


(3.18) 


Since  the  perturbed  solution  is  nearby  the  unperturbed  solution,  we  have 
A=cu^-(o,  |A|  «u)^,  u“-u  =11  |T]  |«|u“|.  ( 


(3.21) 


Substituting  from  (3.21)  into  (3.20),  neglecting  products  of  small  quantities 

c 

and  employing  (3.9)-  (3.11)  and  (3.13),  we  obtain 


which  is  identical  with  Eq. (3.16).  However,  it  should  be  emphasized  that 
although  Eq.  (3.16)  does  not  require  the  existence  of  a con^lete  set  of  ortho- 
gonal  eigensolutions,  Eq. (3.12)  does  require  the  existence  of  such  a set. 

p 

One  consequence  of  this  is  that  second  and  higher  order  perturbation  theory 


cannot  be  obtained  from  the  procedure  presented  in  this  paragraph. 


13. 


From  Eq. (3.16)  [or  (3.22)]  for  the  first  perturbation  in  eigenfrequency, 
we  cam  obtain  equations  for  the  first  perturbation  in  phase  velocity  and  wave- 
number  by  following  the  procedure  outlined  near  the  end  of  Sec. 2 of  Ref. 2. 

■flie  equations  are 

e=A/§^,  6 = A/Vjj,  (3.23) 


where 


(3.24) 


and  v„  and  § are  the  unperturbed  phase  velocity  and  wavenumber  respectively 
of  the  Mth  eigensolution. 

As  an  example  of  the  application  of  (3.13),  (3.15)  euid  (3.16)  [or  (3.22)] 
to  a specific  case  consider  a piezoelectric  solid  with  traction  free  surfaces 
and  for  simplicity  a sufficiently  high  dielectric  tensor  that  the  normal  com- 
ponent of  electric  displacement  can  be  taken  to  vanish.  Let  the  solid  be 
subject  to  static  biasing  stresses,  strains  euid/or  electric  fields  applied  in 
such  a way  as  not  to  affect  the  homogeneous  mechanical  or  electrical  boundary 
conditions  determining  the  normal  modes  of  vibration  of  the  solid.  Since  the 
Mth  piezoelectric  eigensolution  perturbed  by  the  bias  is  for  traction  free 
boundary  conditions  and  vanishing  normal  component  of  electric  displacement, 
we  have 

M M 

(3.25) 


and  under  these  circumstances  from  (3.13)  we  find  that  H takes  the  form 

M 


The  qu2mtities  **L^L  surface  integral  in  (3.26)  are  surface 


perturbation  terms,  which  are  to  be  determined  from  the  Mth  eigensolution  due 
to  the  presence  of  the  bias  as  are  the  terms  in  the  volume  integral  in  (3.26). 


i 


14. 


For  traction  free  and  zero  normal  component  of  electric  displacement  boundary 
conditions  in  the  presence  of  the  bias,  we  have 


N, 


M. 


L L 


(3.27) 


which  with  (3.26)  and  the  divergence  theorem  yields 

»M-J  ®'o' 


(3.28) 


which  is  the  form  taken  by  the  perturbation  integral  in  this  special  but 

important  case.  In  (3.28)  and  ^ take  the  values  given  by  the  Mth  ortho- 

N ^ M 

normal  piezoelectric  eigensolution  ^ presence  of  the  bias  and,  conse- 

quently, from  (3.4)  we  have 

^ M . ^ iM 

*liC  * ‘^lCRV^VjR  ®RLC  ,R’ 

= e,„.g"  ..  - e.-f".,  . (3.29) 


" ®LR'^V,R  “ ®LR^,R  ’ 


the  substitution  of  which  in  (3.28)  yields 


Jr''  M M . n ii'i  , 

[c  - g g-  + 2e  -f  g-  - c f f ] dv  . 
^ ^ lCrtv,r  C,!.  RlC  ,R  C,l  lr  ,l  ,r  o 


M 


- 2.M  2.M 


(3.30) 


Thus,  if  the  piezoelectric  eigensolution  and  bias  are  known,  the  perturlsed 
frequency  u)  can  be  determined  from  (3.15),  (3.16)  and  (3.30). 

4.  perturbations  from  Electroelastic  Solutions 

In  this  section  we  obtain  the  equations  for  perturbations  from  solutions 
of  the  linear  electroelastic  equations  for  small  fields  superposed  on  a bias 

from  a Green*  s function  formulation  of  the  equations  of  linear  electroelasticity. 

* 

TO  this  end  we  first  write  Eqs. (2.9)  and  (2.10)  in  the  form 


J 


15. 


~ o** 

" "v 

" fLVMC'^CjM  '*’ 


(4.1) 


(4.2) 


•^L  “ 2LMC“C,M  '*■ 


and  note  that  the  definitions  in  (2.11)  - (2.13)  hold.  In  (4.1)  and  (4.2)  we 
have  taken  the  liberty  of  introducing  the  extra  quantities  cuid  which 
denote  mechcuiical  amd  electrical  perturbation  terms^  respectively,  and  can, 
in  particular,  represent  the  influence  of  small  material  viscosity.  Prom 
(2.11)j^,  (2.12)^  ^ emd  (2.13)  we  note  that 


iLNMC  ” iMCly’  2^  ~ 2^ 


(4.4) 


Before  presenting  the  Green' s function  formulation  of  the  linear  electro- 
elastic equations  for  small  fields  superposed  on  a bias,  we  show  that  the 
vibrational  eigensolutions  satisfy  an  orthogonality  condition  because  we  need 
this  result  to  obtain  the  full  perturbation  theory.  To  this  end  consider  two 
eigensolutions  of  (4.1)  and  (4.2)  minus  the  perturbation  terms  containing 
and  , one  solution  with  eigenfrequency  (u  and  the  other  with  u)  . Itiey 

L V 


satisfy  the  respective  equations 


V 


IC  + p 0)  u = 0 . 

^V,L  ^ V Y ' L,I< 


(4.5) 


(4.6) 


Prom  (4.5)j^  and  (4.6  )j^  form  the  equation 
I ,~v 


J ' '^v,  L * L * f “ 


(4.7) 


16. 


Performing  the  usual  operations  , eii^>loying  (4.3),  (4.5)2  and  (4.6)2  and  the 
divergence  theorem  in  the  usual  meumer,  we  obtain 


It  should  be  noted  that  (4.8)  cem  readily  Ije  obtained  even  if  intermediate 
surfaces  of  discontinuity  exist.  Clearly,  from  (4.8),  for  homogeneous  boxmdary 


conditions,  we  have 


f P°uf* dv  “ 6 

J ^ Y V o "oi)Vv’ 


(4.9) 


which  is  the  orthogonality  condition  for  linear  electroelastic  vibrations 
superposed  on  a bi2is. 

since  and  are  perturbation  terms,  the  equations  for  the  Green' s 
9 

tensor  for  the  linear  electroelastic  equations  for  small  fields  superposed 


on  a bias  may  be  written  in  the  form 


\y  ” fLYMC®C,M‘^2*®‘Y^,M  ’ 
\ ~ 2LMC°C,  M ^LM^,  M » 


(4.10) 


(4.11) 


(4.12) 


where  equivalent  quantities  are  defined  as  in  the  wording  following  Eq.  (3.7). 
In  (4.10)  - (4.12)  we  have  assumed  that  all  variables  have  a time  dependence  e^' 
We  now  make  the  seune  assumptions  in  (4.1)-  (4.3)  and  in  the  usual  manner, 
from  (4.1)  and  (4.10)  we  form 


17. 


where  h = 1-3  and  all  variables  in  (4.13)  have  spatial  dependence  only.  Per- 

7 

forming  the  usual  operations  , employing  (4.2)^  (4.3),  (4.11)  and  (4.12)  and 

the  divergence  theorem  in  the  usual  manner,  we  obtain 


(4.14) 


Similarly,  from  (4.2)  and  (4.11)  we  form 


(4.15) 


and  performing  the  usual  operations  , utilizing  the  divergence  theorem  and 
employing  (4.1),  (4.3),  (4.10)  and  (4.12),  we  obtain 


?<!■)  - j "^o®’ 

V*  ^ ' 


4 ~4 


(4.16) 


Equations  (4.14)  and  (4.16)  constitute  the  Green's  function  (or  tensor)  formu- 
lation of  the  linear  electroelastic  equations  for  small  fields  superposed  on 
a bias,  it  turns  out  that  for  our  purposes,  although  we  envisage  use  for  (4.14) 
we  do  not  envisage  use  for  (4.16)  because  of  the  particular  type  of  perturba- 
tion problem  in  which  we  are  interested. 

We  now  assume  that  a complete  set  of  eigensolutions  u!|J',  cp  exists  and 


define  orthonormal  eigensolutions  to  our  eigenvibration  problem  by 


(4.17) 


We  now  expand  the  mechemical  displacement  Green' s tensor  G^  and  electric 
potential  Green*  s vector  f in  the  forms 


18. 


= I.  , f = E , 


(4.18) 


where  ^ and  constitute  orthonoxmal  solution  functions  satisfying  the 
appropriate  homogeneous  form  of  (4,1)  - (4.3)  svibject  to  the  appropriate  homo- 
geneous boundary  conditions.  Substituting  from  (4.18)  into  (4.10),  employing 
(4.12)  and  the  homogeneous  form  of  (4.10)  for  every  {i,  contracting  with 


integrating  over  and  utilizing  (4.9),  we  obtain 


Substituting  from  (4,18)  and  (4.19)  into  (4.14),  we  obtain 

where 

"l.  ■ I 

s 

o ^ 

«»  * " 


~ 1LYMC^,M  ^MLY^M' 


(4.19) 


(4.20) 


(4.21) 


(4.22) 


The  perturbation  procedure  is  obtained  from  (4.20)  by  following  the  procedure 
outlined  in  Sec. 3,  after  Eg,  (3.11),  and  the  resulting  equation  for  the  first 
perturbation  of  the  eigenvalue  teJces  the  form 


■ V%  • 


(4.23) 


Again,  as  in  Sec. 3 and  by  following  a procedure  directly  analogous  to  the  one 
employed  in  obtaining  (3.22),  it  may  readily  be  shown  that  (4.23)  can  be  ob- 
tained without  the  use  of  a Green*  s tensor  or  a con^lete  set  of  orthogonal 
eigensolutions. 


19. 
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